Magnon and Hole Excitations in the Two-Dimensional Half-filled Hubbard Model 
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Spin and hole excitation spectra and spectral weights are calculated for the half-filled Hubbard 
model, as a function of t/U. We find that the high energy spin spectra are sensitive to charge 
fiuctuations. The energy difference A(7r, 0) — A(7r/2, 7r/2), which is negative for the Heisenberg 
model, changes sign at a fairly small t/U « 0.053(5). The hole bandwidth is proportional to J, 
and considerably larger than in the t-J models. It has a minimum at [■k/2,%/2) and a very weak 
dispersion along the antiferromagnetic zone boundary. A good fit to the measured spin spectra in 
La2Cu04 at T = IQK is obtained with the parameter values U = S.leV, t = 0.35eV. 

PACS numbers: 75.10.Jm 



I. INTRODUCTION 

Underdoped phases of high temperature superconduct- 
ing materials and the nature of the metal insulator tran- 
sition upon doping a Mott-insulating antiferromagnet 
remain central topics of research in condensed matter 
physics. Some puzzles extend all the way to the un- 
doped stoichiometric insulating materials. Results such 
as the antiferromagnetic zone-boundary magnon excita- 
tions probed in inelastic neutron scattering^'^, the two- 
magnon excitations probed in Raman scattering^ and the 
one-hole excitations probed in angle-resolved photoemis- 
sion spectroscopy^ contiune to surprise us. The question 
of whether some of these anomalies are connected to the 
pseudogap phase of the weakly doped materials remains 
a topic of debate. 

An important question is the extent to which conven- 
tional approaches, based on ordered antiferromagnetic 
phases, can explain the observed spectra and spectral 
weights and to what extent the interpretation of data 
necessitates the introduction of novel ideas such as spin- 
liquids and spin-charge separation. The low energy long- 
wavelength spin excitations of the antiferromagnet are 
well described by the non- linear sigma model^. However, 
the high energy zone-boundary spin excitations neces- 
sarily require a microscopic lattice model. The anti- 
ferromagnetic insulator, without charge fluctuations, is 
represented by the Heisenberg model, and the excitation 
spectrum of this model has been the subject of several 
controlled numerical studies^'^'^. It is clear that the anti- 
ferromagnetic zone-boundary spectrum of La2Cu04 does 
not agree with that of the Heisenberg model. In partic- 
ular, in the Heisenberg model, the magnon energy dif- 
ference A(7r, 0) — A(7r/2, 7r/2) is negative but it is found 
to be positive for La2Cu04. This result will be worse 
if second neighbor antiferromagnetic interactions are in- 
cluded. It has been suggested that one way to reconcile 
the difference is by invoking ring-exchange terms^'^'^'^°, 
which arise due to charge fluctuations^^. 



Here, we present systematic numerical calculations of 
the magnon and hole spectra and spectral weights of 
the Hubbard model as a function of t/U. First, we 
focus on the magnons. Earlier the magnon spectra 
were studied by mean-field theory^^ and by a Quan- 
tum Monte Carlo Simulation combined with the Single 
Mode Approximation^^, neither of which are expected to 
be quantitatively accurate for small t/U . Our calcula- 
tions show that the zone-boundary magnon energies are 
very sensitive to charge fluctuations and the difference 
A(7r, 0) — A(7r/2,7r/2) changes sign at a relatively small 
t/U value of 0.053(5). The magnon spectra of La2Cu04 
and the spectral weights are well described by the Hub- 
bard model as discussed below. 

The calculated hole spectra, on the other hand, are 
qualitatively similar to previous theoretical studies of 
Hubbard and t-J models^**' The hole-bandwidth is 
suppressed at large J7 by a factor of t/U , although we 
find it to still be much larger than in the corresponding 
t-J models. The minimum is at (7r/2,7r/2) with a weak 
dispersion along the antiferromagnetic zone boundary. 
Hence, these results cannot be used to fit the observed 
ARPES spectra in the undoped cuprate materials^. Al- 
though same-sublattice hopping terms can allow better 
fits to the dispersion, the anomalous spectral weights re- 
main harder to explain^. We note that a complete un- 
derstanding of the ARPES experiments may require a 
multi-band model, as well as inclusion of dielectric and 
charging effects. 

To carry out an Ising type expansion^^'^^ for this sys- 
tem at T = we consider the Hubbard- Ising model with 
the following Hamiltonian: 

H = Ho + XHi 

Ho = J/4^(afa| + l) + ^[C/(niT-i)(n^-i) 

(ij) i 
-f/l(-l)Vf] 

Hi = -J2[J{a!a^ + l)/A + t{clci, + h.c.)] 
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FIG. 1: (Color online) Magnon dispersion curve along se- 
lected directions in the Brillouin zone for various values of 
U/t, expressed in units of Jeff = At^ /U. Results for the 
Heisenberg antiferromagnet are shown as a red curve. 




FIG. 2: Magnon energies at two wavevectors along the anti- 
ferromagnetic zone boundary as a function of t/U . 



-/i^(-l)Vf 



where af — — nij, and A is the expansion parameter. 
The Ising interaction J is, in principle, an adjustable pa- 
rameter but here is chosen to be At^ /U . The strength 
of the staggered field h can be varied to improve conver- 
gence. Note that the full Hubbard model is recovered at 
A = 1, at which point the extra terms cancel between Hq 
and Hi. On the other hand for A < 1, there is an Ising- 
like anisotropy in the system, which favors a Neel state 
and induces a gap in the spectrum. The limit A = cor- 
responds to the Ising model, with the usual Neel states 
being the two unperturbed ground states. 



We have extended the linked cluster method"'^^ to the 
spectra of the Hubbard-Ising models. At A = 0, the 
model has a very simple excitation spectrum. Above the 
two ground states, all single spin- flip states (or all states 
with a single hole for the hole spectra) are degenerate 
with each other. We construct an orthogonality transfor- 
mation, which order by order in powers of A, decouples 
these single spin-flip states from the rest of the Hilbert 
space. This procedure amounts to finding the combina- 
tion of one-flip states with others that remains an eigen- 
state as A changes from zero. For a translationally in- 
variant system the resulting block diagonal Hamiltonian 
in the one particle subspace is diagonalized by Fourier 
transformation. 

To set the normalizations for our calculations, we begin 
with the definition for the dynamic structure factor 

1 r°° 

'^"''(q,^) = ^^y^^exp[z(c.t + q.(r,-r,))] 



(1) 



We define the single magnon contribution to the dynam- 
ical transverse structure factor as 

5^^(q,c^)-f 5^^(q,c^) = As5{uj-eM))+B,{ci,uj) (2) 

Here, es(q) gives the dispersion for the magnons and 
As(q) defines the weights for the magnons. The quantity 
Bs{(\,u}) defines the multiparticle contributions. Series 
expansions have been calculated for the spectra ts (q) and 
the spectral weights ^^(q) up to order A^^. 

Similarly, for the hole-excitations, we define the spec- 
tral function 

^(q, ^) = 1^ Yl j exp[i(u;i + q • (r^ - Vj))] 

{c]jt)c.Amt (3) 

We define the single hole contribution to the spectral 
function as 



A{q, Lo) = AhS{LU - e,,(q)) + Bhiq, w) 



(4) 



Series expansions are calculated for the hole spectra e/i(q) 
and the spectral weights Ah{q) up to order A^^. 

The gap in the spectrum closes at A = 1, when spin 
rotational symmetry is restored in the model. This 
causes power- law singularities in certain properties of the 
model^^. Hence, the A = 1 limit needs to be dealt with 
by series extrapolation methods. We use the method 
of integrated differential approximants, well known from 
the study of critical phenomena^^, to calculate various 
properties at A = 1. 

In Figure 1, we show the calculated magnon dispersion, 
in units of JeS = At'^/U, along selected directions in the 
Brillouin zone, for several values of t/U. The results^ 
for the Heisenberg model are also shown. While the long 
wavelength spin-wave velocity is gradually reduced with 
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FIG. 3: (Color online) Fits to the spectra of the material 
La2Cu04: experimental spectra at lOK (open symbols) fit by 
U/t = 8.75±.7, Jeff = 162±3meV (green curve); experimental 
spectra at 295K (solid symbols) fit with U /t = 10.5 ± .7, 
Jeff = 152 ± SmeV (red curve) . 
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FIG. 4: (Color online) Comparison of magnon spectral weight 
g^As for U/t — 10.5 with neutron scattering intensity^ for 
La2Cu04 at 295K. 



increasing t/U, other dramatic changes arise along the 
antiferromagnetic zone boundary. The magnon energy 
at wavevector (tt, 0) at first rises briefly before it begins 
to decrease with increasing t/U. On the other hand, 
the magnon energy at (7r/2,7r/2) decreases sharply with 
increasing t/U . Both of these are plotted in Figure 2, 
where one can see that they cross at a relatively small 
value of t/U = 0.053(5). 

In Figures 3 and 4 we show fits to the magnon spec- 
tra and spectral weight of La2Cu04. Since we have not 
done any finite temperature calculations, we try to fit 
the spectra at different temperatures by effective U and 
t values. We find that the spectrum at lOK is fitted well 



by U/t = 8.75 ±0.7 and Jos = 162±3 meV, while that at 
295K is fitted best by f//i 10.5 ± 0.7 and Jeff = 152 ± 3 
meV. These results suggest that as the temperature is in- 
creased the effective exchange constant decreases whereas 
the effective U/t ratio increases. Assuming that the lOK 
data are essentially at T = 0, we obtain bare parameters 
of f/ = 3.1eV and t = 0.35eV. 

The spectral weights are only measured at 295K, hence 
we show a fit to the calculated spectral weights at the 
larger U/t ratio. In any case, the spectral weights are 
not very sensitive to the U/t ratio. The relative fit 
is excellent. To get a measure of the absolute fit, we 
first note that the integrated one-magnon spectral weight 
over the entire zone was^'' found from experiment to 
be 0.36 ± 0.09, in the normalization where total trans- 
verse spectral weight is 0.5. With this normalization, 
our more accurate calculations for the Heisenberg model 
give an integrated transverse one-magnon spectral weight 
of 0.419(2). The one-magnon spectral weights decrease 
slightly with decreasing U/t, and are also much less ac- 
curate, giving 0.40(8) for U/t — 10. Hence, the results 
agree with experiments well within the uncertainties. We 
should note, however, that this agreement is very much 
dominated by the \/q dependent behavior near the anti- 
ferromagnetic wavevector (where q is the deviation from 
the antiferromagnetic wavevector), which does not de- 
pend much on the microscopic model. The best place 
to look for multimagnon excitations and a more sensi- 
tive comparison with microscopic models is the spectral 
weight along the antiferromagnetic zone boundary. For 
the Heisenberg model, the multiparticle spctral weight is 
largest at (tt, 0), where it is about 40 percent of the total 
transverse spectral weight. We hope our work may moti- 
vate more accurate measurements of multiparticle spec- 
tral weights along the antiferromagnetic zone boundary. 

The data on two other systems of square-lattice anti- 
ferromagnets (CuDC00)2 • 4D2O (CFTD)^! and Cu(II) 
spins of Sr2Cu304Cl2^^ have much larger U/t ratios and 
are well fitted by the Heisenberg model. Allowing U to 
vary, (CuDC00)2 • 4D2O (CFTD) can be well fitted by 
the Hubbard model with U/t = 50 and U = 3.9eV. On 
the other hand, the experimental data on Sr2Cu304Cl2 
have substantial uncertainties and anisotropics, so one 
can not get a reliable estimate for U/t. If one assumes 
the U value for it is about 4eV, one estimates U/t k, 40. 

In Figure 5, we show the single hole excitation spec- 
tra along selected directions in momentum space at dif- 
ferent values of t/U . Plots are also shown for the t-J 
model^'^, with comparable t/ J values. It is evident that 
the band-width scales with J but is much larger than in 
the t-J model. The reason for this is the effective same- 
sublattice hopping generated in the Hubbard model in 
the order t^ /U which is not included in the t-J model. 
In all cases, the minimum of the hole energy remains at 
(7r/2,7r/2) and the dispersion along the line (7r/2,7r/2) 
to (tt, 0) remains weak. As mentioned above, this is in 
contrast to the observed single hole dispersion** in the 
material Sr2Cu02Cl2. Furthermore, our calculated spec- 



4 



'^'---M^ iJ/t=0.25 



U/t=16, J/t = 0.25 




(0,0) (;i,0) (7r/8,Ti/2) (0,0) 



FIG. 5: (Color online) Single hole dispersion in the square- 
lattice half-filled Hubbard model along selected directions in 
momentum space for several values of t/U (squares). Also 
shown as crosses is the dispersion for the t-J model with com- 
parable t/Jeff ratios. 



tral weights at (7r/2,7r/2) and (7r,0) are very similar. It 
has been argued that even small same-sublattice hopping 
terms can significantly change the shape of the disper- 
sion curves and bring them closer to those observed in 
ARPES measurements^, they may also help reconcile the 
measured spectral weights in the undoped cuprates^^. 

In conclusion, we find that allowing for charge fluc- 
tuations by making t/U finite allows us to understand 
the antiferromagnctic zone-boundary excitations in the 
material La2Cu04 very well. However, it appears that 
the simple one-band Hubbard model considered here is 
unable to explain the single hole spectra in the undoped 
cuprate materials. 
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